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In this paper, we study the conductance of the graphene nanoribbons(GNRs) in the presence of 
the Stone-Wales(S-W) reconstruction, using the transfer matrix method. The ribbon is connected 
with semi-infinite quantum wires as the leads. The S-W reconstruction occurs on the edges and 
the interfaces between the electrodes and ribbon. When the reconstruction occurs on the edges, 
the conductance is suppressed considerably if the gate voltage Vg takes intermediate values around 
IFgl ~ to (to being the hopping amplitude of grahene) in both positive and negative energy regions. In 
contrast, if Vg is close to the Dirac point or the band edges, the conductance is relatively insensitive 
to the edge reconstruction. The effect of edge reconstruction become less important with increasing 
ribbon width as expected. The S-W reconstruction occurs also possibly at the interfaces. In this 
case, the reconstruction suppresses identically the conductance in the entire range of Vg for armchair 
GNRs. For the zigzag GNRs, the conductance is strongly suppressed in the negative energy region, 
however the change of the conductance is relatively small in the positive energy region. We also 
analyze the transmission coefficients as functions of the channel index(the transverse momentum 
ky of the leads) for the neutral armchair GNRs with interface defects. Interestingly, there are two 
transmission peaks appearing at ky = tt/S and ky = 2tv/S due to the unit cell doubling. 


1. INTRODUCTION 

Graphene is one of the carbon allotropes with extraor¬ 
dinary electronic properties [Il-Q, that make it prominent 
for potential applications in nano-electronic devices. In 
practice, one need to take the edge effects into account, 
since they could change the band structure and have im¬ 
portant implications for the transport properties of nano¬ 
size graphene. The conventional edges like armchair and 
zigzag have been well studied, however they may not 
be robust against edge reconstructions in some realistic 
circumtancesQ. Indeed, the reconstructed edges have 
already been observed in experiments 

One of the reconstructions is caused by the Stone- 
Wales(S-W) mechanism, which is obtained by 90° rota¬ 
tion of the carbon-carbon bond without changing total 
atom number. In Fig. [T] (a) and (c), we show two typ¬ 
ical continuous edge reconstructions by the S-W mecha¬ 
nism, ac-757 in armchair graphene nanoribbons(AGNR) 
and zz-57 in zigzag graphene nanoribbons(ZGNR), re¬ 
spectively. Without hydrogen passivation, a single ac- 
757 defect contains two adjacent heptagons and one pen¬ 
tagon, while the zz-57 reconstruction consists of alternat¬ 
ing heptagon and pentagon on the edges [ 3 , 0 . In con¬ 
trast to the graphene nanoribbons(GNRs) with the ideal 
edges, the more realistic S-W edge reconstruction breaks 


electron-hole(e-h) sjunmetry and changes the band struc¬ 
tures substantially S 0 , of which the effects on the 
transport of GNRs have attracted attension recently[§- 
P. It was found in Ref. dl that in the wide and short 
GNRs the conductance change due to the S-W edge re¬ 
construction is inconspicuous when the gate voltage is 
close to the Dirac point. While, in another Ref. 0 it is 
shown that in the narrow and long GNRs, the S-W de¬ 
fects at the edges modify the band structure and result 
in strong backward scattering, that suppresses the con¬ 
ductance. These two seemingly different senarios may be 
attributed to the ribbon size, and also the distribution 
of the S-W defects. In this paper we study the trans¬ 
port through GNRs with the S-W defects periodically 
distributed on the ribbon edges. We also study the dif¬ 
ferent effects of the S-W reconstructions when they occur 
at different positions, i.e, the ribbon edges or the hetero¬ 
junctions between the ribbon and the electrodes. 

The paper is organized as following. Our model and 
method are briefly introduced in Sec. 2, and more details 
can be find in Refs.[T^.[T3|. The effect of S-W reconstruc¬ 
tion is given in Sec.3, where we only consider the ac-757 
and zz-57 reconstructions. The conclusion and summary 
are given in Sec.4. 


2. MODEL AND METHOD 


The GNR is put in connection with semi-infinite quan- 
*Electronic address: bugubird'zhang@hotmail.com tum wires as the electrodes which are simulated by the 

lElectronic address: xiaoqunwang@ruc.edu.cn Square lattices. Since the edge reconstruction is in- 
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(a) AGNRs with defects at the edges (b) AGNRs with defects at the interfaces 



FIG. 1: (Color online.) The schematic illustration of S-W reconstruction on the ribbon edges (a, c), and at the interfaces 
(b, d) between the ribbon and the electrodes. The square lattices with width M are the leads, and the current flows in the 
a;-direction. L and Na denotes the length and width(in terms of carbon atom number) for the GNRs, respectively. 


evitable in the fabrication and transfer of GNRs, it can 
occur finally on both the ribbon edge and the interfaces 
between the GNRs and electrodes. Figs. [T] shows the 
geometry of GNRs with ac-757 and zz-57 reconstruction 
on the edges or interfaces. The Hamiltonian defined on 
this network can be expressed in a tight binding form 

" = - E WiTidr-itEGG. (1) 

{ij/j') d 

where Vg is the gate voltage applied on the GNRs and 
Cjj {Cij) is the creation (annihilation) operator at lat¬ 
tice site i,j. j' is the hopping amplitude between 
two nearest neighbors ij and i j . The hopping coeffi¬ 
cients on the carbon hexagon edges is to which is usually 
equal to 2.7eV, while those on the deformed bonds be¬ 
longing to heptagons or pentagons are taken from Refs. 
[l^ as shown in Fig. [1] and Tab. ID For simplicity, the 
hopping parameters in the leads and the GNR-electrode 
interfaces are assumed to be tp. The Fermi surface is 
set at zero energy, and the spin indices of electrons are 
omitted simply for convenience. 

In the leads with width M, the quantum states can 
be labeled by the quantized transverse momenta ky^n = 
mT/{M -b 1) which correspond to different conducting 
channels. According to the Landauer-Biittiker formula, 
the conductance is the summation over all the transmis¬ 




t2 

to 

t4 

t5 

te 

tj 

is 

ti/to 

0.982 

1.189 

0.991 

1.018 

0.982 

0.957 

1.009 

1.189 


TABLE I: The hopping amplitudes for the deformed bonds in 
unit of to, which are taken from Refs. [l^ . ti’s are defined in 

Fig.m 

sion coefficients from channel n to n' 

2e^ ^ 

G = — ( 2 ) 

n,n =1 

r„ Ti.' can be calculated using the transfer matrix method 

[1III3. 

We assume all the S-W defects align in a periodic pat¬ 
tern on the boundaries of the GNRs with width Na and 
length L(see Fig. [Ija, c)), which breaks the translational 
symmetry and the unit cell is doubled. The randomly 
distributed defects are not discussed in this article. The 
spectra of GNRs with periodic ac-757 and zz-57 edge de¬ 
fects are shown in Fig. [2] which are consistent with those 
in Refs. [i,®- It turns out the e-h symmetry is com¬ 
pletely lost due to the existence of 5 and 7 atom rings 
on the edges. As a consequence, for the AGNRs with ac- 
757 defects, the conduction and valence bands touch at a 
positive energy A > 0 for mod{Na,3) = l(see Fig.^a)). 
There appear four new subbands corresponding to edge 
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(a) Defected AGNRs,Na=79 (b) Defected ZGNRs,Na=80 
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FIG. 2: (Color online.) (a) is the band structure of AGNRs 
with ac-757 edge reconstruction, and (b) is that of ZGNRs 
with zz-57 edge reconstruction. Note that the unit cell is 
doubled by the reconstruction, and we take the length of the 
original unit cell in the pristine GNRs as unit, then the mo¬ 
mentum kx belongs to the interval [0,7r). 

states [l^, indicated by colored lines in Fig. [5] (a). For 
ZGNRs in the presence of the zz-57 reconstruction, the 
midgap band is no longer flat, but bends towards the 
valence band(see Fig. [2](b)). Furthermore, there occur 
additional midgap subbands with higher energy. These 
new features of the energy spectra due to the edge recon¬ 
struction might affect the transport properties finally. 

3. TRANSPORT THROUGH AGNRS AND 
ZGNRS WITH S-W DEFECTS 



FIG. 3: (Color online.) The normalized conductance G/M as 
functions of the gate voltage Vg for AGNRs with S-W recon¬ 
struction on the ribbon edges (a) and (b) and at the interfaces 
(c) and (d). The black curves are the conductance Go of the 
pristine AGNRs for comparison. Inset of (b) shows the ratio 
Go/Go as a function of the width at Vg = 0 and ±to- 

We explore in details the transport properties of GNRs 
with ac-757 and zz-57 reconstructions in this section. 
Figures 13] (a) and (b) show the normalized conductances 
Gd /M of AGNRs with periodic S-W defects on the rib¬ 


bon edges as functions of the gate voltage Vg. The over¬ 
all e-h asymmetry of the conductance data is mainly due 
to the odd-number-atom rings at the interfaces between 
the leads and the AGNRs which makes the lattice non- 
bipartitefl^. The conductance is suppressed consider¬ 
ably in the intermediate energy region 0.7to < Vg < 2to 
and —2.7to < Vg < —O.lto by the ac-757 edge reconstruc¬ 
tion, where the trigonal warping effect is important which 
is incompatible with the symmetry of heptagon and pen¬ 
tagon. For Vg close to zero energy or the band edges, Gd 
is almost the same as that of pristine AGNRs denoted by 
Go, as shown in Fig. 131(a), except that the position of the 
minimal conductance shifts slightly from zero towards a 
negative value. This is due to the breaking e-h symme¬ 
try by the ac-757 edge reconstruction, which results in 
the shift of the touch point of the valence and conduc¬ 
tion bands where the D.O.S. as well as the conductance 
is minimal. 

In Figs. |3Ka) and (b), we also show the size effect on 
the conductance. When the ribbon length increases, the 
electrons experience more scattering from the S-W edge 
defects leading to the further suppression of the conduc¬ 
tance as shown in Fig. |3Ka). While, the increase of the 
ribbon width has two effects, one is yielding more con¬ 
ducting channels leading to the linear dependence of the 
conductance on the width, the other is weakening the 
edge effect so that the conductance converges to that of 
pristine AGNRs when the ribbon width is increased, as 
shown in Fig. 131(b) and its inset. 

When the zz-57 defects locate at the interfaces between 
the leads and the AGNR, the heptagons and pentagons 
at the interfaces strongly scatter electrons and suppress 
the conductance in the entire region of the gate voltages 
as shown Figs. m c) and (d). In this case, we take a rela¬ 
tively large ribbon width Na > 200, so that we can focus 
on the reconstruction at the interfaces and ignore those 
on the ribbon edges, i.e., we are liberated to choose the 
conventional armchair edges. In this setup, the conduc¬ 
tance data are obviously insensitive to the size as shown 
in Figs.|3Kc) and (d). 

Next, we consider the transport properties of ZGNRs 
with the S-W reconstruction. Without edge reconstruc¬ 
tion, the lattice of the ZGNR connected with the leads is 
still bipartite unlike the case of AGNRs [l^, therefore the 
conductance data remains e-h symmetric as shown with 
the black curves in Fig. When the zz-57 reconstruc¬ 
tion is presented on the edges, the e-h symmetry is bro¬ 
ken due to the occurrence of heptagons and pentagons, 
and the midgap band is not flat anymore(see Fig.|3Kb)). 
This is reflected in the conductance data near Vg = 0 
in the inset of Fig. |4Ka), where the peak position shifts 
from zero(black curve) toward a positive value(red curve) 
and is enhanced due to the dispersion of the midgap 
band induced by the edge reconstruction. Except for the 
shift of the midgap peak, the zz-57 edge reconstruction 
barely changes the conductance data for small gate volt¬ 
age. This situation holds roughly in the energy interval 
|Vg| < 0.5to and \Vg\ > 2.3to- For the intermediate en- 
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FIG. 4: (Color online.) The normalized conductance G/M as 
functions of the gate voltage Vg for ZGNRs with S-W recon¬ 
struction on the ribbon edges (a) and (b) and at the interfaces 
(c) and (d). The black curves are the conductance Go of the 
pristine ZGNRs for comparison. Inset of (a) shows the data of 
Go (black line) and Gd (red line) with L = 40 and Na = 80 
in the energy region [— 0.2fo,0.2to]. 



FIG. 5: (Color online.) The transmission coefficients in differ¬ 
ent fcy-channels for AGNRs without any reconstructions (a, 
b), with the S-W reconstruction on the edges (c) and at the 
interfaces (d). The vertical dashed lines mark the position of 
27r/3. 


ergy region OMq < \Vg\ < 2.3to, the influence of the edge 
reconstruction is obvious. Figure HKb) shows the width 
dependence of the conductance data, which also implies 
the edge effect is not important in wide ribbons. 

Figures SKc) and (d) show the conductance data with 
the ac-757 reconstruction at the interfaces between the 
leads and the ZGNR. A relatively large width Na = 200 
is taken so that the ribbon edges do not have essential ef¬ 
fect on the bulk transport, and can be choosen as zigzag 
pattern for the simplicity of computation. The e-h sym¬ 
metry is obviously broken. The conductance in the neg¬ 
ative energy region is strongly suppressed by the scatter¬ 
ing from the heptagons and pentagons at the interfaces. 
In the positive energy region, the influence of the ac- 
757 deformation is not strong, and besides the reduction 
of the conductance as expected, there is a small region 
0 < Ig < 0.2to where the conductance is actually en¬ 
hanced. Similar to the case of AGNRs, the conductance 
is insensitive to the ribbon sizes when the S-W recon¬ 
struction occurs at the interfaces. 

Finally, we show the transmission coefficients at Vg = 0 
through different channels ky^n, which is related to the 
states localized at the interfaces of the electrodes and the 
GNRs [T6l - [l8l| . For pristine AGNRs connected with the 
leads, it is the zigzag pattern that occurs at the inter¬ 
faces. As well known there are zero energy states local¬ 
ized at the zigzag boundaries with momenta ranging from 
27r/3 to 47r/3 [l^[23|, which are also involved in electron 
transmission as long as their localization length is com¬ 
parable with the ribbon length. This leads to the trans¬ 
mission peak around ky = 2tt/3 as shown in Figs. [SKa) 
and (b). The peak position is not sensitive to the ribbon 
size, however the intensity is strongly dependent on the 
ribbon length (see Fig. EKa)) since increasing the ribbon 


length reduces the number of the midg ap states extend¬ 
ing from one interface to the other |16l |. 

When the ac-757 reconstruction occurs at the edges 
of AGNRs, the e-h symmetry is broken and the touch 
point of the conduction and valence bands shifts upward 
as shown in Fig. [2] (a). As a consequence, the maxi¬ 
mal transmission peak moves slightly from ky = 27r/3 
toward a smaller value, accompanied by additional satel¬ 
lite peaks as shown in Fig. [5jc). These satellite peaks 
arise because of the boundary condition of the zigzag in¬ 
terface is modified by the zz-57 edge reconstruction. 



FIG. 6: (Color online.) (a) is the transmission coefficients 
in different fcg-channels for ANGRs with S-W reconstruction 
at the interfaces with various lengths L . (b) is the subband 
structures relavant for the transport at zero gate voltage in 
ZGNRs with the S-W edge reconstruction(see also Fig.[2](b)). 
The black dashed lines mark the momenta tt/S and 27r/3. 

When the S-W reconstruction occur at the interfaces, 
there is an additional zero energy transmission peak at 
ky = 7r/3 with smaller intensity, besides the original 
one at ky = 27r/3, as shown in Fig. [S] (d). The posi¬ 
tions of these two peaks deviate from tt/S and 27r/3 for 
small ribbon length L, though they converge to tt/3 and 
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27r/3 eventually as L increases as shown in Fig. |6l[a). 
As expected, this variation coincide with that of the 
zero energy momenta with the width Na of the corre¬ 
sponding ZGNRs with zz-57 edge reconstruction shown 
in Fig. |6l[b). However these momenta given in Fig. [6] don 
not have a naive one-to-one correspondence, because the 
momentum in Fig. El^b) belongs to the folded Bril- 
louin zone(BZ), while ky in Fig. IHlJa) corresponds to the 
unfolded BZ. In fact, for ZGNRs with zz-57 edge recon¬ 
struction, the translational symmetry is broken with unit 
cell doubled, then the two states with momenta k and 
7r-|- fc in the unfolded BZ are mixed leading to new states 
with momentum k in the folded BZ. Therefore in Fig. [6] 
(b), the zero energy state at kx = 27r/3 in the folded 
BZ corresponds to the mixing of two momenta 27r/3 and 
57r/3 in the unfolded BZ, which match the two conduct¬ 
ing channels with ky 2 tt/3 and 7r/3 in the leads, re¬ 
spectively. The same holds for the zero energy state at 
kx = Tr/d in the folded BZ, which corresponds to the mix¬ 
ing of the two momenta tt/S and 47r/3 in the unfolded BZ 
also corresponding to the two conducting channels in the 
leads also with ky ~ 7r/3 and 27r/3. Thus, we obtain two 
transmission peaks at 7r/3 and 27r/3 in Fig. [5] (a). 

4. CONCLUSION AND SUMMARY 

In the presence of edge reconstructions, the band struc¬ 
ture of GNRs change substantially. For example, the e-h 
symmetry may be lost completely and the midgap bands 
may appear in both AGNRs of ZGNRs with their edges 
reconstructed by the S-W mechanism. This could even¬ 
tually affect the transport through GNRs, in particular 


those with small size, which is crucial for designing the 
GNR-based nano-devices. 

When the S-W reconstruction occurs on the ribbon 
edges, the conductance in the intermediate energy 
region around Vg ^ tg is suppressed considerably, in con¬ 
trast, Gd changes little for Vg close to zero energy or 
band edges. When the S-W reconstruction occurs at the 
interfaces between the GNRs and the electrodes, the con¬ 
ductance is suppressed in the entire region of Vg for AG¬ 
NRs, while the conductance of the ZGNRs shows strong 
e-h asymmetry, i.e., it is strongly suppressed in the nega¬ 
tive energy region, but change only a little in the positive 
energy region. Therefore, to pursue high conductance of 
GNRs, the defects should be avoided in heterojunctions 
between the electrodes and the GNRs. 

We also study the transmission coefficients in different 
channels of AGNRs with possible S-W reconstructions 
when = 0. If the zigzag pattern occurs at the inter¬ 
faces, there is only one transmission peak at ky ^ 27r/3, 
with possible satellite peaks in the presence of ac-757 de¬ 
fects on the ribbon edges. However in the presence of 
zz-57 reconstruction at the interfaces, two transmission 
peaks appears at 7r/3 and 27r/3 attributed to the unit 
cell doubled by the alternating heptagon and pentagon. 
Our study may be useful in designing the GNR-based 
nano-devices. 
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